NORDITA-94/6 P 

On primordial magnetic fields* 

Kari Enqvist^ 
Nordita 

Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
February 1, 2008 



Abstract 

A primordial magnetic field could be responsible for the observed magnetic fields 
of the galaxies. One possibility is that such a primordial field is generated at 
the electroweak phase transition because of the fluctuations in the Higgs field 
gradients. I describe a statistical averaging procedure which gives rise to a field 
of a correct magnitude. Another possibility, where the Yang-Mills vacuum itself 
is ferrromagnetic, is also discussed. 
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1 Introduction 



The very early universe is in notoriously short supply of observables that would have 
survided until today It is likely that the baryon number of the universe and the density 
perturbations are one. Another, but a much more speculative possibility, might be the 
magnetic fields of spiral galaxies, the origin of which still largely remains a puzzle. 
The nearby galaxies have magnetic fields of the order of B ~ 10~ 6 G [|TJ, which can be 
deduced from observations of the syncrotron radiation put out by electrons travelling 
through the fields, assuming equipartition of magnetic and particle energies^. Recently, 
a field of a similar magnitude has been observed also in a object with z=0.395 [0. 

The model for galactic magnetic fields most studied is the galactic dynamo |4j], 
where differential rotation and turbulence of the ionized gas amplifies a weak seed 
field by several orders of magnitude. Not much is known about the seed field. As the 
dynamo growth time of the magnetic field cannot be smaller than the galactic rotation 
period r ~ 3 x 10 8 yrs, this gives a lower limit of -B see d^l0~ 19 G on a comoving scale 
of the protogalaxy, about 100 kpc. In the Milky Way and the Andromeda Nebula 
the dynamos appear to be rather weak and the growth time as long as r ~ 10 9 yrs 
which would imply that S see d^-B exp(— t/r) ~ 10~ 10 G. Moreover, in the Milky Way 
the magnetic field changes its direction by about 180° between the Sagittarius and 
the Orion spiral arms ||, and it has been argued || that such a reversal implies a 
stringent lower bound of -B S eed<^10~ 7 G on the seed field. As such a reversal has only 
been observed in the Milky Way, it might not be a generic feature. 

One interesting possibility is that the seed field is truly primordial, with an origin 
that predates nucleosynthesis. In that case the protogalaxy collapsed with a frozen- 
in magnetic field, which enhanced the primordial cosmological field by a factor of 
10 4 Thus at the scale of 100 kpc the dynamo mechanism requires a primordial 
field somewhere in the ballpark of of 10~ 18 G, with an uncertainty of a few orders of 
magnitude. 

From a theoretical point of view, however, the generation of a sufficiently large 
persistent magnetic field in the early universe is rather difficult. There are various 
attempts, relying on phase transitions such as the cosmic inflation or the QCD phase 
transition ||, but the field often comes out to be too small to be of cosmological 
interest. It has been suggested that a large field might actually be generated at the 
electroweak phase transition because of random fluctuations in the Higgs field [|]. If 
one assumes a stochastic, uncorrelated distribution of the Higgs field gradients, or of 
the magnetic field itself, one finds [IT]] today at 100 kpc a root-mean-square field of 
the order of 10 -19 G, which could well serve as the origin of the seed field. This positive 
result is based on calculating the statistical averages along an arbitrary curve. This is 
not the only possibility, but averaging over areas or volumes would produce a field far 
too small to be of relevance for the dynamo effect. 

1 Equipartition may not to be valid for certain irregular galaxies ||. 



In Yang-Mills theories there is also the possibility JTT| that the vacuum is an ana- 
log of the ferromagnet with a non-zero background magnetic field. This is a non- 
perturbative effect, and the resulting field is typically very small. If one is willing, 
however, to go up all the way to the GUT scale one finds that a typical GUT phase 
transition could have given rise to a background field large enough to serve as the seed 
field IO. 



2 EW magnetic fields 

Electromagnetism first occurs when the standard electroweak SU(2) ® U(l)y theory 
is broken down to U(l) em . It is therefore particularly attractive that Vachaspati [[| 
has explained the origin of a primordial field in terms of the cosmological boundary 
condition that all physical quantities should be uncorrelated over distances greater 
than the horizon distance. Since we start with the group SU(2) eg) U(l)y before the 
electroweak phase transition, the resulting electromagnetic field can be constructed in 
a way which is different from the usual F^ u = d^A v — d v A^. The result is f| 

Fa = -i{v!v s -ViV}) , 

v > = Gil/?** ■ (1) 

where <ft is the Higgs field. At the electroweak phase transition the correlation length 
in the broken phase is ~ 1/mw (assuming that the Higgs mass is comparable to mw)- 
The field is thus constant over a distance ~ 1/mw, but it varies in a random way 
over larger distances. Its variation is due to the fact that the field makes a random 
walk on the vacuum manifold of <fi. The problem then is to estimate the field F^ 
over a length scale ~ N/rnw- If iV = 1, it then follows that on dimensional grounds 
Fij ~ rn^y ~ 10 24 G, with probably an uncertainty of ±1 in the exponent. For iV 
large, one should use an approriate statistical argument. The issue is, which is the 
appropriate way to average over the random fields. 

Let us consider random fields walking around in space in a certain number of steps. 
Thus we replace the continuum by a lattice, where the points are denoted by greek 
letters a, . . .. I wish to estimate the magnetic field over a linear scale (which at most 
is equal to the horizon scale). Thus, let us consider a curve consisting of N steps in 
the lattice and define the mean value 

— 1 N 

5=t^E^ , (2) 

i=l 

where B is a component of the magnetic field, and where the lattice points are on 
the curve. 

Now this curve is arbitrary, and we could take any other curve. Let us therefore 
define the average (...), which averages over curves spanning an iV 3 lattice, i.e. over 



all space. Then, for example, 



1 N 



(B) = ^(J2 Bai ) > (3) 
JV i=i 



which means that for each curve with A" steps the mean value B is computed, and 
this is done for a set of curves which span an iV 3 -lattice, and the average is then 
computed. Therefore (B) depends in general on N, but for simplicity of notation we 
shall leave out the explicit reference to this dependence. It should be emphasized that 
the ensemble average (El) takes into account the field value at each lattice point, so 



that the average is really over the whole lattice volume |10 
Similarly, one can define higher moments such as 

_ 1 N 

(B 2 ) = T72 E (B^B*) , (4) 
JV i,j=l 

together with quantities like ((B — (B)) 2 ) . Note that in ([|) the sum is over curves of 
length iV steps of the non-local quantity (B at B aj ). 

In the stochastic variable was taken to be the Higgs field itself which varies 
over the vacuum manifold. Vachaspati argued that the gradients are of order 1/yN, 
since makes a random walk on the vacuum manifold with A0 ~ \/N, and since 
Ax ~ N. Thus Vi is, in a root mean square sense, of the order 1/yN, and hence 
Fij is of order 1/N. Taking further into account that the flux in a co-moving circular 
contour is constant, the field must decrease like l/a(t) 2 , where a(t) is the scale factor. 
Using the fact that in the early universe a goes like the inverse temperature, the field 
was then estimated to behave like (F^)? ~ T 2 /N. For a scale of 100 kpc this leads 
to (Fij) now ~ 10~ 30 G, which is far too small to explain the galactic fields (unless 
there exists some large scale amplification mechanism). One should also point out 
that this assumption presumes that the total magnetic flux through a given surface is 
a stochastic sum of the fluxes through the unit cells of that surface, or in other words, 
that the fluxes through two adjoining unit cells are uncorrelated. Whether this is true 
or not is an open question. 

It is however natural to assume that also the gradient vectors Vi are stochastic, as 
was done in This is because they directly specify whether there is a magnetic 



field or not, whereas this is only true indirectly for the Higgs field itself. Also, the 
vectors Vi are relevant for questions of alignment between neighbouring domains. 

3 Random Higgs gradients 

Consider the expression ([!]) for the magnetic field in terms of the Higgs gradients V,. 
It is convenient to split these fields into real and imaginary parts, 



K(x) = Ri(x) + 27, (x) 



(5) 



where Ri and Ii are real vectors. Let us consider the system at a fixed time. The 
cosmological boundary condition is then that Ri and Ii are random fields. Let us 
make the following assumptions: 

(i) The random fields have a Gaussian distribution. Thus, the mean value of some 
quantity Q is given by 

(Q) = II / ^^~Q e-^W-^)) 3 -^?-^)) 8 , (6) 

where D is a normalization factor defined such that (1) = 1, and A is a measure 
of the inverse width. The quantities Ri and Ii are the mean values of Ri and 
Ii defined along a curve of length N steps. Q Thus, eq. (0) is relevant for a 
3-dimensional world which is an N 3 lattice. 

(ii) I assume that the mean values are isotropic, i.e. = (R 2 ) = (R3) and 
(Ii) = </ 2 > = (Is). 

Assumption (i) is certainly the simplest way of implementing lack of correlation 
of the gradient vectors over distances compatible with the horizon scale, whereas as- 
sumption (ii) is natural as there is no reason to expect any preferred direction. 

It should be noted that the distribution (H) factorizes into an /-part and an i?-part. 
Thus, for any expectation value consisting of J's and -R's one has factorization, 

(Rl 1 ■ -K 1 n--- 1 t) = (X 1 ■ ■ -KlCln ■ ■ -3£> • (7) 

This property turns out to be very useful in computing the higher moments. 

The expectation value of a component Bi of the magnetic field can now easily be 
found. From the expression (|I]) one finds that 

Bi = -eijkFj k = —i EijkVjVk = 2eij k RjI k . (8) 

Thus 

1 JV I N 

Si = ^E5f=2^E« • (9) 

iV 1=1 iV 8=1 

Hence 

2 N 

iV i=i 
2 N 

= ^m(E (RT - (RiWT - (T*» + N(Ri)(h)) ■ (10) 
1=1 

Now, due to the factorization (|7|), the first term on the right-hand side of the last Eq. 
( [TDD vanishes^, and hence 

(Bj) = 26^(1,) (I fc ) = e jlk ((Ri)(l k ) - (Rk)(li)) =0 (11) 

2 Note that this implies that the mean value in a point is assumed to be equal to the mean value 
computed along all curves of length N. Thus the mean values can depend on TV. 
3 Because — (Bi)) = (1%* — (Ik)) = for symmetry reasons. 



because of the isotropy assumption (ii). Consequently the mean value of the magnetic 
field vanishes. 

The second order moment is given by 

(B 2 ) = —Y,(R a R p ■ l a i p - R a i p ■ i a R p ) 

= -^E{(rf)(/^}-(^)(/f/;)} , (12) 

a/3 

using the factorization (^). Now 

(R?R*) = 11/ ^ R?R^~^ 

= Yx 5i ^ P + n/ [W^ + - W<^>] ^ A(iJ7 ~ <7!>)2 (is) 

and similarly for Further we have e.g. 

d?IV „fl _wp-Y_/"5\\2 t-t /" d?B? 



n y ^-i^e-^-^y = n y ^-^-m^-^ +m 

= (Rj) , (14) 



i.e., the mean value in a given arbitrary point /3 on the lattice is equal to the mean 
value computed over all curves. Using Eqs. flI3|) and ([TJ) in Eq. (fTJ) we get 

<^ 2 > = ^x: + ^ca> 2 + w 2 )) + (w 2 a> 2 - cwa» 2 ) • as) 

The first term is 0(N/N 2 ) = 0(1/ N). The last term, being the square of the mean 
value, actually vanishes because of isotropy. Thus we conclude that the rms value of 
the magnetic field scales like 



'^> = y?(^ + i« 7 > ,+ < fl>, »)~ 0( ^ ) • (16) 

The reason for this slow decrease is the fact that isotropy prevents the mean value 
m entering i 
of short range. 



from entering in (Bj) and (B 2 ), and that the correlations of the gradient vectors are 



4 Consequences of the electroweak magnetic field 

Let us now assume that at the time of the electroweak phase transition, a magnetic 
field with a coherence length £ is generated, with a scaling as given by Eq. fllED- Such 
a field evolves according to usual magnetohydrodynamics 

— — = V x (v x B) — er _1 V x V x B, (17) 



where the conductivity o ~ oo in the early universe. Accordingly, the field is then 
imprinted on the charged plasma^. 

At later times the original coherence length is redshifted by the expansion according 

to 

£(*) = ^o • (18) 
a 

The frozen-in magnetic field is also redshifted by the expansion of the universe. Thus 
at later times at the physical distance scale L = iV£ one finds, 




B ims (t,L) = B ^ = Bol?) mm , (19) 



where T 2 t = 0.301 M P / y g*(To) with g* the effective number of degrees of freedom, 
and U ^ 1.4 x 10 3 (n h 2 )- 2 yrs is the time when the universe becomes matter dom- 
inated; for definiteness, we shall adopt the the value Q h 2 = 0.4, which is the upper 
limit allowed by the age of the universe. 

It is not obvious what the coherence length £ actually is. It is likely that it is 
macroscopic and much larger than the interparticle separation. Let me however assume 
for brevity that Bq ~ 10 24 G and £ ~ 1/T. It is then easy to find from Eq. fll9|) the 
size of the cosmological field today. Taking t = 1.5 X 10 10 yrs and L = 100 kpc 
(corresponding to iV = 1.0 x 10 24 ), we find that today the electroweak magnetic field 
at the scale of intergalactic distances is 

B rms = 4 x 10~ 19 G . (20) 

This seems to be exactly what is required for the numerical dynamo simulations to 
produce the observed galactic magnetic fields of the order 10 -6 G. The inherent un- 
certainties in the estimate ( p0|) are: the value of VLqIi 2 used for computing £*; the time 
at which the magnetic field froze, or T ; the actual value of the field B Q . Therefore 
one should view ( ^0|) as an order-of-magnitude estimate only. 

We should also check what other possible cosmological consequences the existence 
of the random magnetic field, Eq. fll9|), may have. Let us first note that the energy 
density ps in the rms field is very small. In the radiation dominated era we find that 
the energy density within a horizon volume V is 

^hr^-KwS^f- (21) 

The horizon distance is r# = 2t so that ps ~ T 5 /Mp p 7 , and the magnetic field 

contribution to the total energy density is negligible. 

In principle, magnetic fields could modify primordial nucleosynthesis, as discussed 

in ||13|| . However, the electroweak magnetic field at the time of nucleosynthesis and 

at the horizon scale is only -B rms ~ 1500 G which is far too small to give rise to any 

modifications of the Big Bang nucleosynthesis. 

4 There is a possible caveat here: in very large magnetic fields such as considered here the velocity 
of the plasma might depend on the background magnetic field. In the following I will neglect this. 



5 Ferromagnetic universe 



There is also another, more exotic possibility for producing magnetic fields, which is 
based on the observation that, due to quantum fluctuations, the Yang-Mills vacuum is 



unstable in a large enough background magnetic field |1JJ. There are indications from 
lattice calculations that this is a non-perturbative result [Q. Such magnetic field 
fluctuations in the early universe could be sufficient to trigger the phase transition 
to a new, ferromagnet-like ground state with a magnetic field made permanent by 
the charged plasma. In this scenario the primordial field is thus generated as a non- 
perturbative quantum effect. 

The new vacuum results provided the /3-function has a Landau singularity: 

dx 



< oo. (22) 



Then the effective Lagrangian has a minimum away from the perturbative ground 
state Tr F 2 = 0, given by 



-g 2 Tr F 2 u \ min = A 4 , (23) 



1 

T 

where A is the renormalization group invariant scale 

A = ^expf- f -^rV (24) 



where \i is a subtraction point associated with the definition of g. 

The condition for the minimum can be realized in many ways. One of them is a 
constant non-abelian magnetic field Bf = £ijkF°" k with a non-zero component only in 
one direction in the group space, and with a length given by 



gVB a B a = A 2 . (25) 

Consider now SU(N) at the one-loop level. We then have the one-loop, zero 
temperature effective energy for a constant background non-abelian magnetic field 
which in pure SU(N) theory reads [IT| 



y(B) = l;B 2 + ^g 2 B 2 [\n^-^) (26) 




2 96tt 
with a minimum at 

gB min = /i 2 exp (- Y^iJ (27) 

and V min = V(B min ) = — 0.029 (gB min ) 2 . Thus the ground state (the Savvidy vacuum) 
has a non-zero non-abelian magnetic field, the magnitude of which is exponentially 
suppressed relative to the renormalization scale, or the typical momentum scale of the 
system. Thus, for example, for SU(2)l at the electroweak scale the vacuum magnetic 
field would be very small. In the early universe, however, where possibly a grand unified 



symmetry is valid, the exponential suppression is less severe. It is also attenuated by 
the running of the coupling constant. For a set of representative numbers, one might 
consider a (susy) SU(5) model with «gut — 1/25 and Tg UT ~ 10 15 GeV, as in the 
super symmetric Standard Model. This yields B ~ 5 x 10~ 8 yU 2 , which turns out to be 
a magnitude which is relevant for the dynamo mechanism. 

In the early universe the effective energy picks up thermal corrections from fermion- 
ic, gauge boson, and Higgs boson loops. In SU(2) these are obtained by summing the 
Boltzmann factors exp(— f3E n ) for the oscillator modes 



E 2 n = p 2 + 2gB(n + 1) + 2gBS 3 + m 2 (T), 



(2f 



where S3 = ±1/2 (±1) for fermions (vectors bosons). In Eq. (28) I have included 
the thermally induced mass m(T) ~ gT, corresponding to a ring summation of the 
relevant diagrams. Numerically, the effect of the thermal mass turns out to be very 
important. 

The detailed form of the thermal correction depends on the actual model, but we 
may take our cue from the SU(2) one-loop calculation, which for the fermionic and 



scalar cases can be extracted from the real-time QED calculation in |fL5| . The result 
is 

(gB) 2 ^, r°°dx 



sv4 



8V£ 



Air 2 

(gB) 



57T 



Y(-l) l+1 f°° ^e"W [ xcoth(x) - 1] 

U. Jo x 

2 00 

: £ 



1=1 



OO rjnr> 

X 3 



X- 



sinh(x) 

where the normalization is such that the correction vanishes for zero field, and 

gBl 2 



(29) 



Kf{x) 



m 2 x 



AxT 2 + gB 



(30) 



where a = f, b stands for fermions or bosons. 

For vector bosons there is the added complication that there exists a negative, 
unstable mode, which gives rise to an imaginary partf]. At high temperatures the 
instability is absent for fields such that gB < m 2 (T), which is the case we are interested 
in here, so that no regulation of the unstable n = 0, S3 = — 1 mode is needed. Thus 
we find IH2l 



5V£ 



(gB) 2 g 



1=1 



00 dx 



x° 



-K b Ax) 



x- 



cosh(2x) 
sinh(x) 



(31) 



At high temperature, the bosonic contributions are more important than the 
fermionic ones. When B ^ T 2 ~ m 2 (T), we find numerically that 5V£ ~ 0.016 x 
(gB) 2 . This gives rise to a small correction to the magnitude of the field at the min- 
imum as obtained from Eq. (BTj) . We may thus conclude that the Savvidy vacuum 
exists for all T. 



Physically the imaginary part is an indicator that the vacuum also contains vector particles pi 



The transition to this new ferromagnet-like vacuum is triggered by local fluctua- 
tions. Charged particles in the primeval plasma generate current j = V x B. The 
typical interparticle distance is L ~ 1/T and a typical curl goes like 1/L so that 
B ~ jL where j is like charge density with one charge in the volume L 3 . Thus the 
Maxwell equations imply that B ~ 1/L 2 = T 2 , indicating that the creation of the 
Sawidy vacuum can take place locally. A constant non-abelian magnetic field, given 
by Eq. (p7|), is then imprinted on the plasma of particles carrying the relevant charges. 
The Maxwell magnetic field B em is a projection in the space of non-abelian magnetic 
fields, and we take it to be of the size comparable to B in Eq. (^%). 

The magnetic flux remains conserved (recall that the primodial plasma is an ex- 
tremely good conductor), and we may write 

-1 ,.2 / «7r 2 \ ,„l 2r , MfGUT) XJ 



B ( r >^"M-nwJUr 3 * 10 G l m ) ' (32) 

where [i ~ T and a(t) is the scale factor of the universe, and the last figure is for susy 
SU(5). This expression is valid because the energy E of the vacuum is redshifted by 
l/a(t). Now in the minimum E is proportional to VB 2 , where V is the volume. Since 
V is proportional to a 3 , we get B ~ 1/a 2 . Hence the magnetic energy per horizon is 
much less than the radiation energy. 

As time passes, the universe undergoes a number of phase transitions. Each of 
these correspond to new types of ferromagnetic vacua, which in general have decreasing 
field strengths. However, the original GUT vacuum has existed for a time that is long 
enough for the plasma to interact with the vacuum field B given by Eq. (^). This 
interaction does not allow for the GUT flux to decrease once it has been created since 
it has become a feature of the plasma, which conserves the flux in the sense that the 
magnetic lines of force are frozen into the fluid. Even if the original field is suddenly 
removed by the creation of a new vacuum, the magnetic field will survive in a perfect 
conductor (see e.g. pp. 186-189 in ref. |T7| ). 

From Eq. (^) we find that the Maxwell magnetic field at t now ~ 10 10 yr is given 
by B now ~ 3 x 10 42 G(tGUTA*)(£*Anoiij) 4 ^ 3 — 10~ 14 G. Such a magnetic field appears 
comparable to what is needed for the seed field in galactic dynamo models. Note also 
that at nucleosynthesis one obtains B ~ 10 4 G, which is well below the nucleosynthesis 
bound on magnetic fields |TB . 



6 Discussion 



The GUT causal domain Iq has today the size of only about 1 m. Obviously, during 
the course of the evolution of the universe, domains with magnetic fields pointing to 
different directions have come into contact with each other. One might think that this 
results in domain walls. However, here it is important that the magnetic flux lines 
follow the plasma particles and cease to be homogenous. If inside each GUT horizon 



a magnetic field line at a certain time passes through two plasma particles, then this 
is true at any later time. The magnetic field thus "aligns" with the plasma. When 
two horizon bubbles collide, the two plasmas rearrange and become one plasma, and 
the same is true for the magnetic field lines, which are part of the new plasma 
The field "realigns" with the "new" plasma and the root mean square of the magnetic 
field remains of the same order as before. Because there are no domains, no random 
walk factor appears at large distances. 

In this argument it is important that -B rms is much smaller than the square of the 
rms momentum, since otherwise the electrical conductivity would depend on B. This 
condition implies that the radius of curvature of a typical plasma particle is very large 
compared to the mean free path, or that the magnetic energy is much less than the 
kinetic energy of the plasma. This certainly is the case in the ferromagnetic universe 
model. 

The size of the field is determined by the scale at which the ferromagnet vacuum 
is created, and the earlier this happens, the bigger the field. If there is a period 
of cosmic inflation, then the relevant field would be created after reheating. If the 
reheating temperature is comparable to GUT scales, the strength of the magnetic field 
would be given as in fl32|) , with interesting consequences for the formation of galactic 
magnetic fields. In this scenario the primordial seed field is thus a relic from the GUT 
era. 

If the origin of the magnetic field is the electroweak phase transition, the situation 
might be different because at scales 1/T the magnetic energy equals to the energy in 
radiation. In such fields the plasma might be trapped by the field, rather than the 
field being imprinted on the moving plasma. Very likely this would result in a domain 
structure, but in the absence of any true dynamical calculation, the details must 
remain unclear. It nevertheless seems clear that a primordial magnetic field would 
have many intriguing consequences, some of which might actually be observable. It 
would of course be of great interest to detect this relic field directly in the intergalactic 
space. 
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